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Abstract In case of chiral Ua(1) symmetry restoration the mass of the t] boson (the ninth, would-be 
Goldstone boson) is decreased, thus its production cross section is heavily enhanced. The t]' decays (through 
one of its decay channels) into five pions. These pions will not contribute to Bose-Einsten correlations, 
thus the production enhancement changes the strength of two-pion correlation functions at low momentum. 
Preliminary results on Bose-Einstein correlation functions support the mass decrease of the t] boson. In 
this paper we propose a method to select pions originating from j] decays. We investigate the efficiency of 
the proposed kinematical cut in several collision systems and energies with several simulators. We prove 
that our method can be used in all investigeted collision systems. 
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1 Introduction 



\ In relativistic gold-gold collisions of the Relativistic Heavy 
- Ion Collider (RHIC) a strongly interacting quark gluon 
plasma is produced [Ij. The temperature of this matter 
may reach values up to 300-600 MeV f2]. At these very 
high temperatures the degrees of freedom are not hadrons 
but quarks or gluons. It is expected that the broken sym- 
. metrics of QCD may be partially restored in this matter. 
\ In a three-quark QCD there is a [/l(3) x Ur{3) chiral 
. symmetry. According to group theory, C/(3) = SU{3) x 
° t/(l), thus chiral symmetry can be written as SUl{3) x 
SUr{3) X Ua{1) X /7y(l). From this SUl{3) x SUr{3) 
' is the flavor-symmetry, which is spontaneously broken to 
' SUv{3). With this symmetry breaking eight low-mass Gold- 
stone bosons are created, associated with the eight pseu- 
doscalar mesons: three pions, three kaons and the rj. The 
, Ua{^) part is also broken, the ninth light Goldstone-boson 
is however missing [3lll]. This puzzle is solved by the 
Adler-Bell-Jackiw anomaly: the Ua{^) part of chiral sym- 
metry is broken explicitly due to instantons tunneling 
between topologically different QCD vacuum states. The 
Goldstone-boson appearing with this symmetry breaking 
is expected to be massive. It is associated with the r]' me- 
son, which has a mass of 958 MeV, significantly higher 
than that of the other 8 pseudoscalar mesons. 

However, in case of chiral symmetry restoration, if the 
[/^(l) symmetry is partially restored, the mass of the r]^ 
is decreased |3lll] . It is an important aspect of this picture 
that the symmetry has to be still partially restored when 
the T]^ mesons are created, i.e. the quark-hadron transi- 
tion has to happen earlier than the chiral transition. This 
seems to be supported by lattice QCD calculations [5]. 



2 Chiral symmetry restoration and mass 
modification of the rj^ 

As mentioned above, chiral symmetry may be partially 
restored in hot QCD matter, and the mass of r]' meson 
might be lower than it's original mass, 958 MeV. How- 
ever, the production cross sections of the light mesons are 
exponentially suppressed by their mass. Hence without 
mass modification roughly two orders of magnitute less r]^ 
mesons are produced than pions. In contrast, decreased 
mass T]' mesons will be created more abundantly. The en- 
hancement of T]' production, according to the Hagedorn 
formula, may be given as [Ij: 



\mn' 



l-d/2 



(1) 



if the mass of the r]' rur^' is decreased to m*, . Here Tcond 
is the temperature of the medium when the rj' mesons are 
created, while d is the effective dimension of the expansion. 
Thus the number of r]' mesons is closely related to their 
mass. 

The thermalized quark gluon plasma is created at RHIC 
in gold-gold collisions [l] roughly 1 fm/c after the collision. 
The matter expands for a time estimated to be around 6- 
10 fm/c and cools down to the quark-hadron transition 
temperature range 150-170 MeV [5j. The hadrons created 
at this point are then propagating freely to the detectors. 
Some of them however decay throughout their journey, as 
is the rj' meson. It has a mean lifetime of 1000 fm/c, thus 
during its life the medium dissolves and the symmetry is 
broken again, so the rj' mass increases again at the expense 
of its momentum: 
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where the starred quantities refer to the in-medium prop- 
erties, while the others to the vacuum properties of the vl . 
Thus the vacuum momentum of the if]' will be significantly 
lower than its original momentum. 

The decay of the vl happens after it regained its vac- 
uum mass. One important decay channel is the decay into 
two leptons, r^' ^ + /~, this is investigated in ref. [6j. It 
turns out that there is an excess in the dilepton spectrum 
at low invariant mass, and this excess might be related to 
the T]' enhancement. There is also a decay mode when the 
T]' goes into an 77 and two pions, and the r\ also decays into 
three pions: 



(7r+ 



- TT 



- TT 



(3) 



and the overall probability of this decay chain is 10% [7J. 
The average momentum of the resulting five pions is 138 
MeV due to the low momentum of the original if]' [8j . Thus 
the T]' decays will have an effect on pion correlation func- 
tions. We will investigate this in the next section. 



3 Two-pion Bose-Einstein correlations 

Final state effects distort two-particle correlation func- 
tions. One of the most important final state effect is that 
of Bose-Einstein correlations. These can be reviewed as 
follows. Definition of the two-particle correlation function 
is: 



C2 (Pl,P2) 



^2 (Pl,P2) 



(4) 



where pi and p2 are the momenta of the two particles, A^i 
and N2 are the one- and two-particle invariant momentum 
distributions. They are defined as: 



A^2bi,P2) = y S{xx,vx)S{x2,V2Wi^2Yd'x2d'xx (6) 

with S{x^p) being the hadronic source function (some- 
times noted as emission function), ^i{x^p) is the one- 
particle wave function and ^1,2(^1,^2,^1,^2) is the two- 
particle wave function. Latter has to be symmetrized in 
case of identical boson pairs. In case of the plain- wave ap- 
proximation (neglecting all other final state interactions), 
we arrive at the following result (see details in ref. [9j): 



^2(^1,^2) 
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S{q,pi)S{q,p2) 
^(0,pi)^(0,p2)' 



(7) 



where q = pi — P2 and S{q^p) is the Fourier transformed 
of the source function (the Fourier transformation is in 
X ^ q). Introducing K = {pi -\-p2)/2 and taking into 
account pi — P2, we get: 



C2{q,K) = l 



\S{q,K)\^ 
\S{0,K)\' 



(8) 



In the core-halo model [9j, the hadronic source is di- 
vided into two parts: a core and a halo. The core consists of 
the primordial particles and decay products of very short 
lifetime resonances. This part of the source has a small 
size: roughly 10 fm at maximum. The halo consists then 
of decay products of long lived resonances, such as 77, rj' 
or Kg. The halo hadrons are created very far from the 
core (note the very large lifetime of the previously men- 
tioned particles). When measuring correlation functions 
however, due to finite momentum resolution of the detec- 
tors very small momentum differences cannot be resolved, 
i.e. pairs with such similar momenta are regarded as one 
by the detectors. In the Fourier transformation, large sizes 
correspond to small momenta, thus the halo correlations 
are not seen in measurements. Therefore we introduce the 
core source function, and denote it by Sc{q-,K) (and use 
the H subscript for the halo part). It holds that 



S{q,K) = Sc{q,K) + SH{q,K) 



(9) 



However, for measureable q values (at least several MeV 
values), S{q^K) = Sc{q^K) (again, large size of Sh cor- 
responds to small width of Sh)- Let then the number of 
particles in the core be A^c, number of particles in the 
halo be Nh- Clearly, for g = 0: 



Nc + Nh 
Nc 



SciO,K) 



(10) 



thus finally 



C2{q,K) = l + \ 



\Sc{q,K)\^ 
\Sc{0,K)\' 



(11) 



(5) with A being 



Nc 



Nc + Nh 



Thus 



C2{q^0,K) = l + K 



(12) 



(13) 



The A parameter is thus the q ^ extrapolated value of 
the correlation function C2 , and it depends on the ratio of 
the core to the halo. Thus if the mass of the r]^ is decreased, 
more of it are produced (see in the previous section), their 
decay pions will be also be enhanced in number, so the 
halo will be larger. This means, that the A parameter is 
decreased [8j. Hence r]' mass and A value are connected. 

It was found [HITOlITT] that the A parameter is indeed 
decreasing at the kinematical domain of r]' decay pions. 
However, it is not experimentally proven that the rj' decay 
pions are causing the decrease. In this paper we investigate 
a method to kinematically filter out pions from rj^ decays. 
If applied to the experimental sample, in case of an r]' 
mass modification the A decrese will vanish. 
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4 Kinematical domain of pions from ml 
decays 

Our method is based on the invariant mass of pions from 
the decay chain of eq. ([3]). The invariant mass of pion pairs 
in this decay is: 



mL = (^i+^2) -(Pi+P2f 



(14) 



with E\ and E2 being the energy of the pions, pi and p2 
their three-momentum. Using E'^ = ^ m? we get 



=2ml + 2^m2 -^pj^Jml^pl- 2pip2 coscp (15) 

where (f is the angle between the two pions. If being in 
the rest system of the Er^f = m^/ holds, and due to 
momentum conservation = —pi —p2 is also true, so we 
get 



+ H + P2 + 2pip2 cos (16) 



which can be substituted into eq. ([T5]) . The 77 can take 
most of the energy if pi = p2 (thus (p = 0), and it has 
the least energy if (p = tt. This yields a lower and an up- 
per bound for mf^^ , and the result for the interval will be 
0.078-0.168 GeV^/c^. Similarly for the second part of the 
decay chain the same calculation can be done, and we get 
0.078-0.166 GeV^/c^. This can be checked in our simula- 
tions, and the intervals could been verified. Based on the 
simulations, we chose the 0.075-0.171 GeV^/c^ interval for 
all pairs. Thus this interval provides a selection method of 
T]' decay pions. However, it is not the most effective, since 
a significant fraction of all other pion pair is also in this in- 
terval. Thus we also checked the invariant mass of all four 
pions from this decay. It falls in the 0.43-0.69 GeV^/c^ 
interval, which, together with the two-pion invariant mass 
cut, yields an effective method of kinematical selection of 
T]' decay pions. See example plots on fig. (TJ from 200 GeV 
center-of-mass energy p+p collisions, simulated with HI- 
JING 1.411. 



5 Method of tagging r]' decay pions 

Our method, based on the previous section, is the follow- 
ing. Let us have for example a 7r+ from the experimen- 
tal sample. Let us then select any tt", 7r+,7r~ to form a 
quadruplet. We then check if the invariant mass of the +,— 
pairs is in the above mentioned interval, and also check if 
the invariant mass of the quadruplet is in the mentioned 
interval for quadruplets. If for our starting 7r+ there exist 
three other pions (one same sign and two opposite sign) 
with which all invariant mass criteria are fulfilled, we tag 
this pion "found", as quite probably it comes from an r]' 
decay. A different method is, when we start from a 7r+,7r+ 
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Figure 1. Invariant mass distributions for y/s =200 GeV p+p 
collisions from Hijing simulations. The first (second) plot shows 
the rriinv distribution of all pairs (quadruplets). The last three 
plots show the rriinv distribution for pairs and quadruplets 
coming from r]' (or from the 77 of the same decay). 



pair and select an opposite sign pair to form a quadruplet. 
Then we check if the quadruplet and the two +,— pairs 
also fulfill the invariant mass criteria, and if yes, we tag 
this pair "found". This latter method is expected to have 
one advantage compared to the single particle method: the 
probability of "finding" a non-77' pair than a non-7^' pion 
is expected to be smaller. If applied to experimental data, 
tagged pairs or pions can be removed from the data thus 
creating a sample poor in r]' decay pions. 

In simulations, we can determine if the pair of the par- 
ticle comes from an r]' decay, so the efficiency of the selec- 
tion method can be tested. Whether using the pair or the 
single particle method, we can form four different groups 
of them: 



a) Comes from an rj' and fulfills the rriinv criteria 

b) Comes from an 77' and does not fulfill the m. 

c) Does not come from an 77' and fulfills the m. 

d) Does not come from an 77' and does not fulfill the m. 
criteria 



mv criteria 
inv criteria 



Here fulfilling the rriinv criteria means that, with our pair 
or particle, a quadruplet can be formed that fulfills all 
ruinv criteria. Ideally all rj' pions go into the first group, 
while all others fall into the last group. Let us call the 
number of pions (or pairs in the other method) in the 
four groups A/'^, A^5, Nc and A^^^, respectively. Then the 
number of all pions (pairs) is Na + + A^c + while all 
pions (pairs) from an rj' are Na-\- N^, and the ratio of the 
number of all pions (pairs) from an 77' over the number 
of ah others is {Na + Ni,)/{Nc + Nd) before the filtering. 
This ratio is N^/Nd after applying our method. Both ratios 
are important, since they are directly connected to the A 
parameter. The "goodness" of our method is basically the 
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Figure 2. Geometrical cuts used in the simulations. PHENIX 
has the smallest acceptance, thus the efficiency is the worse 
there. 



double ratio {Na + Nb)/{Nc + Nd) over N^/Nd, as this 
shows whether the cleaned sample contains less r]' decay 
pions. 

Note that the kinematic acceptance of our detectors 
clearly distorts our method. It is possible that not all four 
pions from an r]' are detected, thus the required quadru- 
plet cannot be formed in the sample. In this case not all 
pions (pairs) will be found (and filtered out), this is what 
we call efficiency (and can be calculated as Na/{Na-\- Ni^). 
Also, if the sample is significantly large, there will be 
a high probability of random quadruplets to fulfill the 
rriinv criteria. This we call loss, and can be calculated as 
Nc/{Nc + Nd)' This causes our experimental sample to 
be smaller, thus we will lose statistics. With proper align- 
ment of the rriinv intervals these effects can be minimized. 
Goal of present paper is however to investigate the effi- 
ciency and loss connected to our the method. A similar 
method was investigated in ref. [12j for e+e~ collisions. 
We test the method in p+p and Au+Au collisions, at sev- 
eral center-of-mass energies. 



6 Results 

We used two simulations to test our method: Pythia 
(version 8.135) and Hijing [H] (version 1.411). In the lat- 
ter, proton-proton and gold-gold collisions could also be 
analyzed, while Pythia was used only in case of proton- 
proton collisions. We also simulated the geometric accep- 
tance of the detectors. In case of the 200 GeV RHIC en- 
ergy, we used the geometry of STAR and PHENIX detec- 
tors, while in case of 14 TeV energy, we used the geometry 
of ALICE and CMS detectors. See details on fig. [H 

In the next part, we will give the following three num- 
bers for each type of simulation: 



Efficiency of finding rj' pions: 
Loss (found non-77' pions): 
Change of rj' fraction: 



Na^m 

Nc 
Nc^Nd 
m / Ng + N, 
Nd N.^Nd 



(17) 
(18) 
(19) 



The optimal value for efficiency is 1 (in this case we could 
cut out all pions coming from rj' mesons), for loss (in 
this case we kept all non-7^' pions) and for the third, the 
purification ratio (if it is zero, then after the cut, there are 
no rj' decay products at all). 

We generated 1 000 000 p+p ^/s = 200 GeV events 
with both Pythia and Hijing. Geometrical cuts largely in- 
fluence our method, but it is working in all cases. However, 
in case of PHENIX cuts, the efficiency is very low, because 
a large fraction of pions are not detected, thus we can't 
"find" them to form quadruplets. See fig. [3] for details. We 
also generated 10 000 p+p ^/s = 14 TeV events. Here due 
to the larger average pion number the efficiency is much 
larger, and the purification ratio is very good. See details 
on fig. Hlfor details. At both energies, the pair method is 
better than the single particle method. We finally gener- 
ated 100 Au+Au sJsnn = 200 events with Hijing. Here 
only the pair method was working, as essentially all single 
particles are found, due to the very large statistics. See 
fig. [5] for details. Note that if we "find" a pair or particle 
that is not coming from an 77', then it will be a loss for 
us, as it decreases our experimental sample. For example 
if the loss is 50%, then the sample is reduced by a factor 
of 2, so the statistical errors will be increased by a factor 
of a/2. 



7 Summary 

We investigated a method of rejecting r]' decay products in 
an experimental sample from high energy collisions. Such 
a method, if used experimentally, would help to validate 
signs of partial chiral symmetry restoration. Our method 
proposes kinematic cuts in the invariant mass spectrum 
of pion pairs and quadruplets. The basic idea is to take 
a given particle or pair, complement it to a quadruplet 
with any other random pions, and check if they fulfill the 
kinematic criteria. If there is such a complementation for 
the given pair or particle, that specific pair or particle can 
be tagged as coming from an rj' . We find that our method 
is working for several systems and several energies. We 
used to simulations for cross-checking purposes. The most 
important system is that of ^/snn =200 GeV Au+Au 
collisions, the pair version of our method is working there. 
This method can thus be used in an experimental analysis 
to search for partial chiral symmetry restoration and the 
modification of the r]' mass. 
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Figure 3. Results from >/s = 200 GeV p+p collisions. The Figure 4. Results from ^/s = 200 GeV p+p collisions. The 
method is working in all cases, but efficiency is very low in method is working in all cases, 
case of PHENIX cuts, for the single particle method. 
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Figure 5. Results from y/s — 200 GeV Au+Au collisions. The 
method is working only in the pair method, as loss is 100% in 
the other case. We do not give the purification ratio here as due 
to our 100% efficiency it always zero (or cannot be calculated 
if loss is also 100%). 



